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ABSTRACT 


Equations  governing  the  coupled  thickness-shear  and  flexui’al 
vibrations  of  isotropic » elaatlc  plates  are  solved  for  the  case  of 
a circular  disk  with  free  edges.  The  results  of  computations  are 
given  for  the  resonant  frequencies  analogous  to  those  which  would  be 
excited  in  an  AT-cut  quartz  disk  with  full  electrode  coating.  A 
marked  difference  is  observed  between  the  frequency  spectrum  of  the 
disk  and  the  spectrum  of  cylindrical  vibrations  of  a rectangular 
plate.  The  difference  is  due  to  the  presence  of  thickness-twist 
modes  of  motion 5 in  the  disk,  vrtiich , along  with  the  thickness-shear 
modes,  are  coupled  to  the  flexural  modes. 


Ijitmiuction 


This  la  the  seventh  in  a series  of  papers  on  vibr-ations  of  plates. 
Hie  first  paper^-  contained  a derivation  of  isotropic  plate  eouations 
which  accoranodate  the  phenomenon  of  coupling  between  flexural  and 
thickness-shear  modes  of  motion.  These  equations  were  later  extended 
to  take  into  account  crystal  plates  , piezoelectric  crystal  plates, 
plates  of  varying  thickness,  plates  with  incomplete  electrodes  and  plates 
with  specially  shaped  electrodes.  (The  first^  and  second^  papers  are 
referred  to  as  I and  II  in  the  sequel.)  The  solutions  of  the  varloius 
equations,  given  as  Illustrative  examples  in  these  papers,  were  all  con- 
fined to  cyllrjdrical  motions,  that  is,  to  one-dimensional  problems. 

We  return,  now,  to  the  simplest  (isotropic)  form  of  the  equations 
to  solve  a two-dimensional  problem.  We  have  in  mind  the  antisymmetric 
modes  of  motion  that  are  excited  in  an  AT-cut  quartz  disk.  However, 
the  additional  matriematical  complexities  intro<iuced  by  the  circular 
shape  of  the  plate  make  it  advisable  to  consider,  first,  the  corres- 
ponding motions  of  an  isotropic  disk.  It  seems  likely  tiiat  the  modes 
of  motion  and  the  frequency  spectrum  of  the  isotropic  disk  will  be  found 
to  differ  insignificantly  from  those  of  the  quartz  disk.  Ihis  is  because 
the  two  controlling  parameters  are  the  shape  of  trie  plate  and  the  ratio 
of  the  plate-flexure  modultis  to  the  thickness-shear  modulus.  Ihe  former 
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is  the  same  for  the  two  oases,  arvl  the  latter  has  been  made  identioal . 
in  the  computations,  simply  hy  choosiiig  Pcicscn’s  ratio,  for  the  iso- 
tropic material,  to  give  the  modulus  ratio  appropriate  to  the  AT-cut  of 
quartz.  The  lateral  elastic  constants,  it  is  true,  are  not  quite  the 
same,  but  the  discrepancy  is  small,  and  its  influence  on  the  motion  is 
probably  slignt* 

We  give,  here,  the  exact  solutions  of  the  equations  for  the  iso- 
tropic , free  disk  and  the  results  of  detailed  computations  of  the  fre- 
quencies in  the  neighborhood  of  the  thickness-shear  frequency  of  an 
infinite  plate.  The  frequency  sjjectrum  of  the  disk  is  found  to  have  a 
different  character  from  that  obtained  in  II  for  the  cylindrical  vibra- 
tions of  a rectangular  plate.  At  first  glance,  it  appears  that  the 
thickness-shear  fundamental  and  its  overtones  are  multi-valued,  i.e,, 
corresponding  to  the  fundamental  and  to  each  overtone  of  thickness-shear 
in  the  rectangi  ‘'ar  plate,  there  appears  to  be  more  than  one  resonant 
frequency  in  the  disk.  Upon  closer  examination,  it  develops  that  the 
additional  resonances  are  due  tc  the  presence  of  a thickness-twist  mode 
and  its  overtones,  in  which  the  displacement  is  a twist  about  an  axis 
normal  to  the  plate.  The  complete  frequency  spectrum  is  the  result  of 
coupling  between  thickness-shear  and  flexure  and  between  thickness-twist 
and  flexure. 


Plate  Euuationa 


In  I and  II , the  plate  uquations  were  deduced  from  the  three- 
dimensional  equations  of  elasticity,  referred  to  a rectangular  coordinate 
system.  Tne  resulting  equations  may  be  transformed  to  a polar  coordinate 
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systoffi,  appropriate  to  the  circular  disk,  but  it  is  simpler  to  rederive 
the  plate  equations  in  the  new  syetom  of  coordinates.  Thus,  the  assumed 
form  of  the  displacements,  in  place  of  l(lO)  or  II (5),  is 


^ %(r,  e,t)  ^ y(r,0,t) 

whero  r and  9 are  polar  coordinates  defined  by  x.  £-  - 

i.e.,  tlie  ^ axis  is  nonn«l  to  the  plane  of  the  plate. 

The  derivation  of  the  plate  eqmtions  proceeds  as  before.  Thtis, 
corresponding  to  1(14,15)  or  Il(lO),  we  have  plate-stress  equations  of 
motion 
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where  jO  is  the  dsnsity,  h is  the  thickness  of  the  plats  and  the  plate- 
stress  components  fxre  given  in  terms  of  the  plate-displacements  by 


^ (i/;.  t-  II ) 


where 


p r £ii/iz  (t~  -y^)  , £ y V are  Young's  vnodulus,  the  shear 

TU  Zn  ^ 

modulus  and  Poisson's  ratio,  respectively,  and  At-  = 17”"/ /Z  ■ 

If  (l ) were  substituted  in  (2 ) the  restilting  plate-displacement 
equations  of  motion  would  be  those  corresponding  to  I(l6)  or  Il(ll).  How- 
ever, in  the  present  case  it  is  more  convenient  to  work  with  the  alterna- 
tive formulation  of  the  eqiAations  given  by  1(62,63).  In  polar  coordinates 
again  omitting  a factor  these  becar.e 
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Solution  of  Equations  of  Motion 


In  view  of  the  inodes  of  motion  of  interest  in  the  present  case, 
the  appropriate  solutions  of  (5)  are 

A,  J,  ( r) 

v - A.  J.(  S^r)  CX3S  B 


= Aj^J^(uir)  -S//7  B 


where  (xy)  is  Bessel’s  function  of  the  first  kind,  of  order  unity,  and 
Ai  f Ai  f ^3  arbitrai’y  constants.  These  are  the  modes  of  circum- 
ferential order  unity.  There  are  also  modes  of  higher  and  lower  circum- 
ferential order  in  v^iich,  for  example,  'yj^  - cos  n6, 

but  we  are  not  concerned  with  them  here. 

For  the  vibrations  of  a plate  with  free  edges , the  boundary  condi- 
tions are 

^ ^ ^ 0 at  r = a (ll) 

where  ci  is  the  radius  of  the  disk.  Using  (l),  (4)  and  (lO),  Equations  (ll) 


become 


[( Z . & ; j;  (-.fa)  - ^ (<r,  - OA, 
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[S,aJjJ;a)  -ZJ,  (<r,a)]  (<r,  - i ) A, 

(12) 

+ [ ( 2 - ^ coa  Z <'coa)J  Aj,  ^0 

Is, a Z(S,a)  - J.(S,a)]cr,  A, 

+ [4a  Jo(JzO)- 

+ J,(o,.a)  A3  = O 


Freouencv  Equation 

Equations  (12)  comprise  a system  of  three  homogeneous,  algebraic 
equations  on  the  unknowns  A,  , . Ihe  vanishing  of  the  determinant 

of  their  coefficients  yields  an  enuations  for  the  determination  of  the  fre- 
quency p through  4 ’ H'ls  secular  equation,  expressed  in  a 

form  convenient  for  desk-calculator  ccanputation , reads; 
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Now^  dT7>^  for  all  values  01  p O while  ^ according 
as  7*  ^ A”  ’ A"  " /A  is  the  frequency  of  the  first 

thickness-shear  irode  of  an  infinite  plate.  Hence,  /2.  vlll  be  real  or 
iinaginary  according  as  7 . Equation  (IIS'!  ^ « a useful  form  for 

p p . In  the  range  p , let  /5  - » so  that  the  frequency 

equation  reads 
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and  I^(i(.)  ^ T,  modified  Bessel '»  functi one  of  the  first  kind. 

Prom  the  relation  /3=^?, 4T  Equation  (6.1,  explicit  expres- 


sions for  the  frequency  may  be  obtained.  Thusj 
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Further  V from  the  relation  ^ ^ a.  , 

djh-^ 

where  c/  is  the  plate  diameter, 

Gopputatiop  of  Freauenciee 

Sqxiations  (13)}  (14)  and  (15)  constitute  the  solution  of  the  prob- 
lem. For  a given  material  the  values  of  V and  are  fixed.  Choosing  e 

0 

value  of  ft  or  ft^  fixes  fpjB  by  (14)  and  determines  an  infinite  set  of 
roots } ^ * of  (13).  For  the  particular  [h  or  each  of  ttiese  roots  deter- 
mines a ratio  } trius  giving  an  Infinite  set  of  al/h  corresponding 

to  every  value  of  ft /p  ■ 

For  the  isotropic  plate}  the  ratio  of  the  thickness-shear  modulus 
to  the  plate  flexure  modulus  is 

^ 7T"  ( I --y)  / 

while}  for  the  .".T-cut  of  quai-tZ} 

- 

^ /z(c„  - c\/c^^) 
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f»ft 

ft<ft 
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where  the  eleetle  constants  are  refen-ed  to  the  coordinates  of  the  (rotated) 
AT  section.  For  this  cut  of  quartz,  ^ = 0.2^3.  Hence,  an  isotropic  plate 
with  V = 0.312  gives  the  same  value  of  j as  an  AT  quartz  plate.  The 
computations  were  made  for  this  value  of  . 

To  facilitate  the  task  of  finding  the  roots  y of  Equations  (13  ), 
for  a series  of  values  of  ^ j tables  were  computed  of  the  functions  ; 

J^(x}/J,(k)  G(. 01)100 

0(. 01)10,  10(. 02)16, 

16(.1)20,  20(1)37 

to  six  decimals.^  The  subsequent  numerical  work,  carried  out  on  a desk 
calculator,  yielded  the  •'esonant  frequencies  of  vibmtion  depicted  in  Fig.  1„ 

The  flerural  modes,  in  the  lower  part  cf  Fig.  1,  are  those  of  even 
order  r = 2, 4?6  ...  , and  are  sinllar  to  those  in  the  lower  part  of  Fig.  1 
of  II.  However,  the  upper  parts  ( f>/i^  >/  ) of  the  two  figures  are  differ- 
ent. For  example,  in  the  case  of  the  rectangular  plate  the  tl'.ickness-shear 
fundamental  couples  with  every  even  flexural  mode  in  turn,  but,  in  the  disk, 
the  thickness-shear  fundamental  appears  to  couple  only  with  alternate,  even, 
flexural  modes.  As  a result,  the  thickness-sliear  fundamental  appears  to 
be  double-valued. 

The  apparent  anomaly  is  resolved  when  it  is  recognized  that  there 
ij  pi'esent  a set  of  thickness-twist  modes  in  addition  to  thickness-shear 
and  flexural  modes.  Mathematically,  the  additional  modes  derive  from  the 
solenoidai,  retarded  potential  H . Their  uncoupled  frequencies  may  be 

5 The  computation  of  the  tables  of  Ja(<)  / was  performed  by  3.  Foley 
on  equipment,  made  available  by  the  Watson  Scientific  Computing  Laboratory 
of  the  International  Business  Machine  Corporation.  Both  tables  are 
stored  in  the  files  of  the  Department  of  Civil  Engineering,  Columbia 
University. 
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obtained,  appro:;ctir;iitelj’’5  by  setting  y-0  , In  the  equations,  so  that  H 
is  the  only  potential  left.  Then  we  have 
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where 


--  ( - I ) J,( cor)  - .1  (ivr) 


Applying  only  the  boundary  condition  - 0 on  r ^ a , ye  obtain  the 

secular  eqimtlon  J^(cL'a)-C>  , or 

The  roots  of  (16),  in  conjunction  with  (14)  and  (15),  give  the  approx- 
imate frequency  spectrum  ( vs.  ) of  the  thickness-twist  modes. 

The  curves  arfe  asymptotic  to  - / and  are  very  nearly  tangent  to  the 

alviost  horizontal  portions  of  the  ciirves  in  Fig.  1.  Ihat  is,  the  thickness- 
twist  inodes  have  frequencies  very  close  to  those  of  the  corresponding 
thickness-shear  modes.  The  effect,  of  two  such  rjolgrit>oring  families,  on 
the  coupling  with  the  flexiiral  modes  is  illustrated  In  Fig.  2.  In  Fig.  2a 
is  shown  the  usual  spectrum  resulting  frcxn  the  coupling  of  a family  of  odd 
thicknesa-sheai'  modes  ( n - 1,3  •••  ) ar^d  even  flsx\xral  modes  ( r = 20,22  ...  ). 
In  Fig.  2b  a family  of  thickness- twist  modes  ( ~ -,3  ...  ),  closely  par- 

alleling 'bnc  thickness -shear  mofiee , is  added.  If  the  result  were  to  be 
interpreted  as  coupling  between  flexure  and  mhickness-shear  alone,  it  would 


appear  that  the  thlckness-shear  mcxies  are  double-valued  and  couple  only 
with  eltei'nace  flextire  modes.  This  would  be  especially  so  if  the  coupling 
were  strong  enough  to  smooth  out  th*'  kinks  in  the  almost  horizontal  por- 
tions of  the  cumres,  as  Is  the  case  in  Fig.  1,  However,  when  the  presence 
of  the  thickness-twist  modes  is  taken  into  account,  there  are  no  anomalies. 

It  should  be  observed  that  the  case  described  here  is  a degenerate 
one  of  coupling  of  three  families  of  modes,  because,  although  there  Is 
coupling  between  thickness-shear  and  flexure  and  between  thickness-twist  and 
f icxui-e , there  is  no  coupling  between  thickness-shear  and  thickness-twist. 

An  example  in  which  there  is  full  coupling  of  the  three  families  will  be 
described  in  a subsequent  paper. 
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Fig.  1:  Resonant  frequencies  of  modes  of  vibration  ol  ciTTiimlerenbiai 

order  unity  of  a circular  disk  ( = 0.312)  as  computed  from  Eqs.  (13  ) > (U)^ 

and  (15).  The  ordinate  is  the  ratio  of  the  resonant  frequency  ( » ) to  the 
frequency  ( p ) of  the  first  thickness-shear  mode  of  an  inflni'^e  plate..  The 
abscissa  is  the  ratio  of  the  diameter  { d ) to  the  thickness  ( h ) of  the  disk. 


Fig.  2:  (a)  Freq>ift>icy  spectrujn  for  couj’lJng  betwean  thickness-shear 

and  flexural  modes. 

(b)  Frequency  spectnjur,  for  coupling  between  thickness-sheer 
and  flexiiral  modes  and  between  thickness-twist  and  flex- 
ural modes,  with  no  coupling  between  shear  and  twist  modes. 


